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I.  Introduction. 

In  lot-by-lot  acceptance  sampling  by  attributes  each  item  of  a 
sample  drawn  from  a lot  of  manufactured  items  is  classified  simply  as 
defective  or  non-defective.  A random  sample  is  drawn  from  the  lot  and 
the  lot  is  either  accepted  or  rejected  depending  solely  upon  the  number 
of  defectives  in  the  sample.  When  the  classification  into  defective 
arid  non-defective  is  made  on  the  basis  of  the  measurement  of  a variable 
quality  characteristic,  and  the  decision  to  accept  or  reject  the  let  is 
a function  of  these  measurements  (as  opposed  to  the  number  of  defectives) 
such  inspection  procedures  are  known  as  sampling  inspection  by  variables. 
Since  inspection  by  variables  makes  greater  use  of  the  information 
concerning  the  lot  than  does  inspection  by  attributes,  whenever  the 
testing  of  the  individual  items  involves  measurement,  variables  plans 
require  smaller  sample  sizes  to  furnish  the  same  degree  of  protection 
than  do  attributes  plans.  Variables  sampling  plans  pertain  to  a single 
quality  characteristic,  and  it  is  visually  assumed  that  measurements  of 
this  quality  characteristic  are  independent,  identically  distributed 
normal  random  variables . Such  an  assumption  will  be  made  throughout  the 
paper. 

Sampling  inspection  by  variables  in  divided  into  three  categories, 
known  standard  deviation  plans,  unknown  standard  deviation  plans,  and 
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average  range  plans.  Known  standard  deviation  plans  are  baaed  upon  the 
sample  mean  and  the  known  standard  deviation.  Unknown  standard  deviation 
plans  are  based  upon  the  sample  mean  and  the  sample  standard  deviation. 
Average  range  plans  are  based  upon  the  sample  mean  arid  the  erage  range 
of  k subgroups  of  r items  (the  sample  size  is  then  kr). 

The  purpose  of  this  paper  is  to  present  a matching  collection  of 
variables  sampling  plans  based  on  known  standard  deviation,  unknown 
standard  deviation,  and  average  range.  Th:oa  plans  are  classified 
according  to  sample  size  and  Acceptable  Quality  l-evel  (AQL).  The  prob- 
ability of  acceptance  at  the  AQL  was  done  to  be  consistent  with 
MIL  STD  105A. 

The  plans  presented  below  have  the  following  properties  which  differ 
from  those  plans  already  published: 

1.  Each  OG  curve  shown  represents  a variables  sampling  plan  cased 
on  known  standard  deviation,  unknown  standard  deviation,  and 
average  range.  In  other  words,  if  the  user  chooses  an  OC  curve, 
he  has  at  his  disposal  the  choice  of  fhe  three  types  of  plans, 
guaranteeing,  for  all  practical  purposes,  the  same  protection. 

2.  The  usual  acceptance-rejection  procedure  in  lot-by-lot  sampling 
inspection  by  attributes  is  to  draw  a random  sample  and  accept 
the  lot  if  the  number  of  defectives  is  sufficiently  small. 

This  i3  equivalent  to  the  procedure  of  estimating  the  percentage 
defective  in  the  lot  by  its  best  estimate,  namely,  the  percentage 
defective  in  the  sample,  anc*  accepting  the  lot  if  this  estimate 
is  small.  Similarly,  sampling  inspection  by  variables  can  be 
viewed  in  the  same  manner.  In  this  paper,  the  acceptance-rejection 
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procedure  is  based  upoD  estimating  the  percentage  defective  in 
the  loti/  with  acceptance  occuring  if  arid  only  if  the  estimate 
is  not  too  large. 

3.  All  existing  sampling  procedures  for  two-sided  specification 
limits  have  the  disadvantage  that  the  probability  of  accepting 
a submitted  lot  with  given  percentage  defective  p does  not 
depend  on  p alone,  but  is  also  a function  of  the  actual  lot 
mean  ^ or  equivalently,  on  the  actual  division  of  the  per- 
centage defective  into  components  lying  above  and  bslow  the 
specifications  U and  & For  this  reason,  it  is  not  possible 
to  compute  a single  OC  curve  giving  the  operating  characteristics 
of  a two-sided  plan.  Such  plans  do  not  yield  a constant  prob- 
ability of  acceptance  for  a given  percentage  defective,  but 
rather  a spectrum  of  probabilities.  Computing  OC  curves  for 
all  possible  divisions  of  p results  in  a band  of  curves.  The  „ 
two-sided  procedure  presented  is  consistent  with  (2).  The 
percentage  below  the  lower  specification  is  estimated  and  the 
percentage  above  the  upper  specification  is  estimated.  The  lot 
is  accepted  if  and  only  if  the  sum  of  these  estimates  is 
sufficiently  small. 

This  procedure  has  the  property  that  the  band  of  OC  curves 
due  to  all  possible  divisions  of  the  percentage  defective  is  so 


1/ 


u 


The  estimates  used  for  the  known  and  unknown  standard  deviation  cases 
are  the  uniformly  minimum  variance  unbiased  estimates  of  the  percentage 
defective  in  the  lot. 

The  two-sided  regions  for  unknown  standard  deviation  plans  in  Sampling 
Inspection  by  Variables . by  Bowker  and  Goode,  McGraw-Hill,  1952,  are 
equivalent  to  those  presented  here.  Both  pieces  of  work  stem  from 
research  conducted  at  the  Applied  Mathematics  and  Statistics  Laboratory, 
Stanford  University,  during  the  summer  of  1951. 
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narrow  as  to  be,  for  practical  purposes,  a single  curve. 
Furthermore,  it  then  follows  that  the  CC  curve  of  the  re- 
lated-one-sided te3t  very  closely  approximates  the  OC  band  of 
the  two-sided  test. 

II.  General  Inspection  Criteria. 

Associated  with  each  inspection  characteristic  are  the  design 
specifications.  . If  only  an  upper  specification  licit  is  given,  the  item 
is  considered  defective  if  its  measurement  exceeds  h.  If  only  a lower 
specification  is  required,  the  item  is  considered  defective  if  its 
measurement  is  smaller  than  L.  If  both  upper  and  lower  specification 
limits  are  specified,  the  item  is  considered  defective  if  its  measurement 
either  exceeds  U or  is  smaller  than  L. 

If  the  percentage  defective  of  a submitted  lot  is  known,  no  sampling 
inspection  is  necessary  to  determine  whether  or  not  the  lot  is  to  be 
accepted.  If  the  percentage  is  sufficiently  small,  the  lot  is  accepted, 
otherwise  rejected.  Since'  such  knowledge  about  the  percent  defective  is 
rare,  a logical  procedure  is  to  estimate  this  percentage  defective  fr  ja 
a sample,  and  accept  or  reject  the  lot  on  the  basis  of  this  estimate. 

A sampling  plan  is  then  described  as  consisting  of  the  sample  size  nj 
a method  of  estimating  the  percentage  defective,  and  a maximum  allowable 
estimated  percentage  defective  allowing  acceptance  of  the  lot,  p*.  If 
an  upper  specification  limit  is  given,  the  percentage  above  this  limit, 
is  estimated  from  the  sample  of  si&e  u.  If  ^ 4p  , the  lot  is 
accepted.  If  a double  specification  limit  is  given,  both  and  are 
computed.  If  p ■ s p*,  the  lot  is  accepted. 

In  the  ensuing  anions  sampling  plans  based  on  known  standard 
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deviation,  unknown  standard  deviation,  and  average  range  will  be  presented. 
These  plans  will  be  indexed  according  to  code  letter^  and  AQL.  The  OC 
curves  for  the  corresponding  plans  (sane  code  letter  and  AQL)  based  on 
known  standard  deviation,  unknown  standard  deviation,  and  average  range 
are  essentially  the  same  and  are  presented  in  Figures  1-16. 


III.  Known  Standard  Deviation  Flans. 

In  this  section,  it  will  be  assumed  that  the  standard  deviation  <P 
of  the  characteristic  measurements  is  known.  A complete  set  of  sampling 
plans  indexed  according  to  AQL  and  code  letter  is  given  in  Table  1. 
Operating  characteristic  curves  for  these  plans  are  provided  in  the 
appropriate  figures. 

Before  presenting  the  general  inspection  procedure  for  known 
standard  deviation  plans,  some  discussion  about  two-sided  specifications 
is  in  order.  For  this  case  information  is  available  about  the  percentage 
defective  in  the  lot  before  a sample  is  drawn.  There  is  a sufficiently 
small  value  of  which  is  a function  only  of  the  AQL  that  assures  in- 
coming quality  worse  than  the  AQL.  Consequently,  before  sampling, 
should  be  calculated  and  the  lot  rejected  (with  no  sample  taken)  if  this 
value  is  smaller  than  the  minimum  value  of  on  the  grounds  that  the 
incoming  quality  is  poorer  than  the  AQL.  Mini  mum  values  of  can  be 
found  in  Table  1.  Intuitively,  this  implies  that  the  combination  of 
narrow  specification  limits  and  large  standard  deviation  insure  against 
the  submittal  of  good  quality. 

The  code  letters  run  from  B to  Q,  omitting  A.  Again  this  was  done 
so  that  there  would  be  somewhat  of  a semblance  between  these  tables 
and  the  corresponding  ones  in  MIL  SIS  105A. 


The  inspection  procedure  is  as  follows: 


..n 

£-xi 
i-1  1 


n 


rU-Xx 


1.  Draw  a random  sample  of  n items  and  compute  x 

2.  a.  For  an  upper  specification  limit  compute  ■ (~~)y 

whore  v * . 

■J  n-l 

C 7 

b.  For  a lower  specification  limit  compute  C£  * (-~-)v 

, f n 

where  v * 1 — r-  . 

^ n-l 

c.  For  a double  specification  limit,  compute  both 

CU  “ (“^)v  and  » (^^)v  if  is  greater  than  the  min. 
Otherwise,  the  lot  is  rejected  before  a sample  is  drawn. 

3.  Enter  Tabic  II  with  and/or  C£  and  read  out  p^  and/or  p^ 

i / 

whichever  is  applicable^' 

A # 

4.  For  an  upper  specification  limit  accept  the  lot  if  py  £ p . 

For  a lower  specification  limit  accept  the  lot  if  ^ p*. 

A t 

For  a double  specification  limit  accept  the  lot  if  pp  ♦ pT  S.  p 

Example . 

The  specified  minimum  yield  point  for  certain  steel  castings  is 
55000  psi.  The  standard  deviation  is  known  to  be  (P  ■ 3000  psi.  A 1# 

AQI-  plan  is  to  be  used  with  a sample  size  of  6.  The  yield  points  of  the 
sample  specimens  are 

62000;  61000;  68500;  59500;  65500;  63900. 

The  following  are  computed  from  the  data: 

1.  x - 63400. 


2b.  C£  - (^^^')(1.0 95)  » 3.07. 
3.  From  Table  II  p^  - .107??. 


^ It  is  shown  in  Section  VI-3  that  these  estimates  are  the  uniformly 
minimum  variance  unbiased  estimates  of  the  true  percentage  defective. 
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4.  From  Table  I p is  2.57^.  Hence  the  lot  is  accepted. 


IV.  Unknown  Standard  Deviation  Plans. 

In  this  section  it  will  be  assumed  that  the  standard  deviation  of 
the  characteristic  measurement  is  unknown.  A complete  set  of  sampling 
plans  indexed  according  to  AQL  and  code  letter  is  given  in  Table  III. 
Operating  Characteristics  Curves  for  these  plans  are  provided  in  the 
appropriate  figures. 


The  inspection  procedure  is  as  follows: 


1.  Draw  a random  sample  of  n items  and  compute  x 


a 

5>i 

i«l  1 


and  s 


y (x.-x)‘ 

i-i  x 


2.  a.  For  an  upper  specif ication  limit,  compute  C * , 

u s 

b.  For  a lower  specification  limit,  compute  C^  - . 

c.  For  a double  specification  limit,  compute  both 

C - &£  and  CT  - . 

Us  L s 

3.  Enter  Table  IV  with  C^  and/or  C^  and  read  out  p^  and/or 

whichever  is  applicable^/ 

4.  For  an  upper  specification  limit  accept  the  lot  if  p*. 

For  a lower  specification  limit  accept  the  lot  if  p,  sp  , 

For  a double  specification  limit  accept  the  lot  if  + ^ S p*. 


Examnle. 


The  specifications  for  electrical  resistance  of  a certain  electrical 

component  is  U - 660  ohms  and  L - 620  ohms.  A .U%  AQL  plan  is  to  be  used 
with  a sample  of  size  10,  The  data  are  as  follows: 


y 


It  is  shewn  in  Section  VI-4  that  theses  estimates  are  the  uniformly 
minimum  variance  unbiased  estimates  of  the  true  percentage  defective. 
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639,  640,  650,  647,  66id;  637.  652,  643,  657,  649. 


The  follcvring  are  computed  from  the  data: 


1.  x ■ 647.6 


647.6-620 

“XoT~ 

660-647.6 

e.ol 


* -j  >15.38 

- 3.44 

- 1.54. 


3.  From  Table  IV 


6.04. 


a 

Pl‘° 

Pjj  - 5.31# 

5 - $L  + Pu  " 5-31*- 


4.  From  Table  III 

p*  - 1.30#. 

Since  p > p , the  lot  is  rejected,  furthermore,  for  possible  review 
board  evidence,  the  best  estimate  of  the  percentage  defective  in  the  lot 
is  5.31#.  Furthermore,  it  is  evident  that  all  the  trouble  lies  in  the 
upper  specification  limit. 


V.  Average  Range  Flans. 

In  this  section  it  will  be  assumed  that  the  standard  deviation  of 
the  characteristic  measurement  is  unknown.  A complete  set  of  sampling 
plans  indexed  according  to  AQ1  and  code  letter  is  given  in  Table  V. 
Operating  characteristic  curves  for  these  plans  are  provided  in  the 
appropriate  figures. 

The  inspection  procedure  i*  as  follows: 

1.  Draw  a random  sample  of  n items  and  group  the  measurements 
into  k subgroups  of  5 (where  n ■ 5k)  with  the  exception  of 
sample  sizes  of  3,  4,  and  7.  For  these  sample  sizes  use  a 
single  subgroup  of  3,  4,  and  7 respectively.  For  each 
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3. 

4. 


subgroup  compute  the  range  (R. ).  Determine  the  average 
k 1 _ n 

range  R » R.  )/k.  Calculate  x » (y~~  x.  )/n. 

i*l  1 i-1  1 

a.  For  an  upper  specification  limit,  compute  Jjj  • (-g-)h4=/ 

bo  For  a lower  specification  limit,  compute  C£  ■ 

c.  For  a double  specification  limit,  compute  both 


C”  - (£p)h  and  C£  - 

Enter  Table  VI  with  Cjj  and/or  C£  and  read  out  py  and/or 
p^  whichever  is  applicable. 

For  an  upper  specification  limit  accept  the  lot  if  ^ < p*. 

For  a lower  specification  limit  accept  the  lot  if  pL  < p . 

/N  A # 

For  a double  specification  limit  accept  the  lot  if  p.j  + p^  < p . 


Example. 

The  maximum  temperature  of  operation  for  a certain  device  is  specified 
as  180°.  A 4*  AQL  plan  is  to  be  used  with  a sample  of  size  15.  The 
sample  items  have  operative  temperatures  of 


178,  175,  174,  158,  172  R][  - 20 

177,  166,  172,  167,  163  \ - 14 

174,  173,  162,  182,  170  R,  - 20  . 

The  following  are  computed  from  the  data: 

1.  x - 170.87. 

2b.  2.394  - 1.21. 

3.  - 11.08* 

4.  From  Table  VI 

p*  » 9.61*. 

Since  p^  > p the  lot  is  rejected. 


2/ 


h is  a factor  found  in  Table  V which  is  a function  of  the  sample 
code  letter. 


_L 


T»bl»  i'or  fat  1m ting  tha  Lot  Parcantaga  Dufacrtlve  far  Plana  Based  oi.  Known  Standard  IXnrUtion 
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VI-1.  General  Theory. 

D.  Blackwell  [1]  has  shown  that  if  x has  density  p^(x),  g is  any 
unbiased  estimate  of  9,  and  T is  a sufficient  statistic  for  £,  then 
E(g|T)  is  also  unbiased,  and  furthermore,  has  a variance  no  greater  than 
that  of  g.  Lehmann  and  Scheffe  [8]  have  extended  this  result  and  have 
shown  that  if  T is  complete^/  then  every  estimable  function^/  h(0) 
possesses  an  unbiased  estimate  with  uniformly  smallest  variance  and 
this  estimate  is  the  unique  unbiased  estimate  of  g($)  which  is  a function 
of  T. 


It  has  already  been  shown  [8]  that  the  sufficient  statistics  for 
the  normal  distribution  are  complete.  Furthermore,  it  is  evident  that 
there  exists  an  unbiased  estimate  of  the  fraction  of  a normal  population 
lying  outside  a fixed  interval,  e.g.,  the  fraction  of  independent  observa- 
tions in  a sample  from  this  normal  population  which  lie  outside  this 
interval.  Consequently,  the  fraction  of  a normal  population  lying  out- 
side a fixed  interval  (p)  is  an  estimable  function  and  has  an  unbiased 
estimate  with  uniformly  smallest  variance,  and  this  estimate  is  the 
unique  unbiased  estimate  of  p which  is  a function  of  the  sufficient 
statistics. 


The  idea  of  basing  the  two-sided  acceptance  procedure  on  the  uniformly 
minimum  variance  unbiased  estimate  of  the  fraction  defective  stems 
from  research  carried  out  at  the  Applied  Math,  and  Stat.  Lab.  by  a group 
including  Professors  A.  H.  Bowker,  K.  Rubin,  and  L.  E.  Moses ; 

2 / 

-1  A statistic  T (more  properly  a family  of  distributions  of  T)  is  said 
to  be  complete  is  E^f(T)  * 0 on  all  S£r  implies  f(T)  ■ 0 except 

possibly  on  a 3et  N for  which  p^(N)  ■ 0 for  all  £6r. 

3/ 

— Given  a random  variable  x with  density  p^(x),  a function  g($)  is  said 
to  be  estimable  if  there  exists  a function  h(x)  such  that  E^hCx)  ■ g($). 
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?fe  are  concerned  with  estimating  the  parameter 

■>  (z- /<)2 

U _ i C 

p ■ 1 - f — ~ — e * (T^  dz  . 

L i2TT^ 

Two  cases  will  be  considered;  namely,  (1)  ^ unknown,  cr>  known  (2)  both 

/ 

and  <r>  unknown. 

Let  xi,x2>"  ■ •>5Si  i^®?®11^^  random  variables  from  the  above 

normal  distribution.  Define  'p‘i  as  the  usual  attribute  estimate  of  the 
fraction  defective  i.e. , the  ratio  of  the  number  of  defective  items  to 
the  sample  size.  It  is  evident  that  '£'•  is  unbiased.  Hence  it  follows 
from  Blackwell  [1]  and  Lehmann  and  Scheffe’s  [8]  results  that  p » ECp* \T) 
is  the  unique  uniformly  minimum  variance  (UMV)unbiased  estimate  of  p, 
where  T are  the  sufficient  statistics  for  the  normal  distribution. 

Since  'p'  i3  the  sum  of  independent  identically  distributed  random  vari- 
ables taking  on  the  values  0 and  1,  it  is  evident  that  p - E(p’\T)  is 
equivalent  to  E(^|T)  where  jTis  defined  as  follows:  Let  y be  any  one 

of  the  observations  (x^,x2>. . . ,xn)  3ay  x^. 

(1)  p(y1,x2,x3,...,xn)  - 0 if  L < y < 0 

• • >xj1)  " j-  otherwise. 

VI-2.  The  Estimate  When  the  population  Variance  is  Known. 

If  the  observations  are  drawn  from  a normal  population  with  unknown 

2 X 

mean  ut  and  known  variance  , then  x - — x,  is  a sufficient 

n i-i  i 

statistic.  In  this  case  p * E(p}x)  is  the  UMV  unbiased  estimate  of  p. 
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If  p is  the  estimate  (l)  then 


p(x)  “ E($|x)  “ Pr-^'p  “ 1)^'}  “ ^ - pr-^L  < 7 < U)x  j ■ 1 - j*  dy 


jl.  n\,x; 


where  g(y,x)  is  the  joint  probability  density  of  y and  x,  and  h(x)  is 
ine  probability  density  of  x. 


n x. 


Consider  the  joint  probability  density  of  y and  x’  * 


■?  mO 


-Uy-  m)2  ♦ (n-iKx'-^;-) 

Jn^T  2<nfc  L ' 1 

2 e 
2TT  a* 


The  transformation 


i - Xs=U£l  ♦ JL 
n n 


leads  to 


g(y^) 


n 

•Jn-1  2"IT 


e 


and  division  by 


h.(x) 


nr 

■■ — 1 — e 

■J  2TV  cf1 


n 

o 

2o>“ 


results  in 
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zLzJlL  , rz.  I e 2(n-l)^ 

h(x)  n_1  r^TTcT' 


Therefore 


P(x) 


L ^ n_1  | 2tt  & 


2(n-l) o 


(7-x)2 


dy 


r n 1^-x 

1 nil  . 

/ 


GO 


-GO 


I2TT 


2 dt  + J e 

ces*1*" 

>/n-l  <7 


2 


dt 


» 1-3  • The  Known- ir<  Acceptance  Cr:.t  crier  - 

The  acceptance  procedure  is  formulated  as  follows:  Accept  a lot 

A * ♦ 

of  items  if  in  the  sample  pip  where  p is  so  chosen  that  if  the 

population  percentage  defective  is  p,  i.e.,  if  the  portion  of  the 

population  lying  outside  (U,L)  is  p,  then  the  probabi.13.ty  of  acceptance 

will  be  Lp.  It  is  shown  that  in  the  one-sided  case  this  is  equivalent 

to  the  well-known  ana  widely  used  procedure:  Accept  if  x < U-k 

Let  K be  defined  by 

e 


GO 

r e 


t r 
2 


|2~TT' 


€ 
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In  the  one-sided  case,  when  L »*  -oo,  then 


cc 


/ . 

n~  U-x 

ln-1 


* 

S P 


we  must  have 


I5M  >K 

-4  n-1  <r> 


rv» 


or 


<U-fS 

— n 


K <P 
P* 


I ^—1 

If  we  take  k * J K ^ the  OG  curves  of  the  two  acceptance  procedures 
will  be  identical. 

For  the  case  where  there  is  a double  specification  limit  and  the 

standard  deviation  is  known,  it  is  evident  that  if  < 2 , the 

2 


incoming  value  of  p,  the  percentage  defective,  is  greater  than  the  AQL. 
Consequently,  the  lot  should  be  rejected  without  a sample  being  drawn. 
The  OC  curve  of  this  two  sided  test  takes  the  form  of  a band,  the  lower 
bound  of  which  is  the  one  sided  OC  curve.  For  incoming  quality  which  is 
no  worse  than  the  AQL  the  upper  bound  corresponds  to  equal  division  of 


the  defective  below  the  lower  specification  limit  and  above  the  upper 
specification  limit.  The  band  here  is  relatively  wide,  but  this  is 
desirable  since  this  gives  "good"  quality  a better  chance  of  acceptance. 
For  quality  worse  than  the  AQL,  the  above  restriction  on  eliminates 
certain  divisions,  including  equal  division,  so  that  the  band  quickly 
tends  to  the  one  sided  curve.  A diagram  of  the  OC  band  is  as  follows: 
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Pci 


VI-4.  The  Estimate  cf  p When  the  Population  Variance  is  Unknown. 

If  the  observations  are  drawn  from  a normal  population  with  unknown 

2 _ 

mean  ^ and  unknown  variance  <r  , then  the  pair  of  sample  values  x,  and 
(x^-x)2  are  sufficient  statistics.  In  this  case  $(x,S2)  - E(^]x,S2) 
where 


n x. 


U A, 

- E — 

f-z  n 


i-1 


n 


=^2 


“ YL  (x. -x) 
i-1  1 


If  p is  the  estimate  (1)  then 


p(x,S2)  - Pr£p  - l|x,S2|  - 1 - 


~ ~2  \ 

jy Li*;0  J. 


h(x,S2) 


dy 


_ 1A  _ 

— A.V 


where  fCy^S*)  is  the  joint  probability  density  of  y,  one  of  the 

observations  in  the  sample,  the  sample  mean  x,  and  the  sample  sum  of 

2 2 ,2 
squares  S*,  and  h(x,S  j is  the  joint  probability  density  of  x and  S . 

It  is  well  known  that  h(x,S“)  is  given  by 


h(SE,S2) 


n-' 


n /- 


r-  ,02x  2 2 o 

wn  (S  ) e 


n-1 

f2T f ^P(~)(2cr2)  2 


o 

s~ 


_ 9 

To  find  f (ypc,S  ) consider  the  Joint  density  of  the  sample.  This 

may  be  expressed  as  the  Joint  density  of  the  mutually  independent  sample 
_ 2 

statistics  y,  x* , and  S'  where 


x» 


n x± 

h** 


i (x^x'  )2 

i»2  1 


(2) 


O' 


\ 


/ 

/ 


n-1 


n-1  2 O' 

■ e 


2\ 

(*'y-)2  V 


ott  rr> 


v *■ 
\ 


1 


/ 


y 


\ 

s»2\ 


) . 


Now  let 


- 1?  - 


x« 


s'2 " s2  - ?£i  (iE-y)2 

7 “ 7 • 


Under  this  transformation  the  expression  (2)  becomes 


(■ a\ 
/ 


n-1  f fry*)2 

2*2l  n"1 


\2i 


where 


K 


n 

n-2 

IS  21T<7'(2CT2)  2 P(~) 


and  the  variables  are  subject  to  the  restrictions 


-oo  < y < oo 
0 S S2  < oo 


x-y 
S 

Dividing  the  expression  (3)  by  h(x,S2)  results  in  the  conditional 
_ 2 

density  of  y,  given  x and  S . 


-IS- 


C'*) 


nr 

>i  n-l 


rt2^) 

jTrr'f"#) 


i {i  - <¥>2 


n-4 

2 


If  new  the  transformation  z - ^ ^ -Tn^I  ^ mac*e  ^he  density  of 

the  random  variable  z is  obtained: 


(5) 


gU) 


P(°-2) 

P(H=2)  p(S=£) 


z 


(|  - D-i  (|  -D-i 

(1-z) 


0 < z < 1 . 


This  will  be  recognized  as  a symmetrical  Beta  distribution  with 
parameters  ^ ~ I.  Hereafter,  g(z)dz  will  be  denoted  by  d^(-^=  - 1). 
Since  “ Pr  -|V  > U(x,S^)  ♦ Pr-^y  < l|x,S^)  then  in  terms 

of  z 


VI -5*  The  Unknown-  Acceptance  Criterion. 

Just  as  in  the  case  where  the  population  variance  is  known  the 
acceptance  procedure  is:  Accept  if  p(x,S  ) < p where  p is  so  chosen 

that  if  the  population  fraction  defective  is  p,  the  probability  of 
acceptance  will  be  L^.  It  is  shown  that  in  the  one-sided  csse;  i,e., 
where  there  is  only  one  specification  limit,  say  U,  this  procedure  is 
equivalent  to  the  well-known  and  widely  used  test  x * ks  < U,  where 


- 19  - 


s ■ S/Jn-1. 


' p*  * 

Let  be  defined  by  j d/5  - 1),  then  p S p if  and  only 


if 


n U-x 


2 >|  n— 1 S -{* p* 


and  since  s * S//n-l  this  is  equivalent  to 


5 4 TT  (1-2/^P*'s  * u • 


If  k is  taken  equal  to  7^  (1-2/4^)  the  procedures  p < p and 

• * ■ / Jr 

x •*  ks  < U will  have  identical  OC  curves.  In  this  caper  the  acceptance 

a U-x  n_-3£ 

criteria  p based  on  the  statistic  — — - are  tabled  rather  than  ■%-,  to 

S o 

conform  with  the  common  use  of  the  square  root  of  the  -unbiased  esstimator 

, 2 
of  O''  . 

For  the  case  where  there  is  a double  specification  limit  one  of  the 
authors  made  numerical  investigations  [11]  of  the  OC  curves  for  various 
uivisiviiS  ox  tne  p or  can  cage  defective.  Xt  was  evident  that  the  band 
was  so  narrow  that  for  all  practiai  purposes  it  can  be  assumed  to  be  a 
single  OC  curve,  i.e.,  the  OC  curve  of  the  one  sided  plan. 

VI-6 . The  Acceptance  Procedure  and  Estimate  Baaed  on  the  Range. 

Since  the  pair  of  sample  statistics  (x,2),  where  2 is  the  sample 
range  or  average  ranged  is  not  sufficient  for  the  normal  distribution 
when  the  mean  and  variance  are  both  unknown,  no  uniformly  minimum 
variance  estimate  of  p which  if  a function  x and  2 can  be  derived.  In 
this  paper  the  following  estimate  of  p(x,S)  is  used: 


In  this  paper,  whenever  pn»«»<vo®  the  subgroup  size  is  taken  as  5> 
so  that  R is  the  average  range  of  m subgroups  of  S (the  sample 
consisting  of  Sm  observations).  For  sample  sizes  of  3,  4,  and  7, 
S is  taken  as  the  range  of  the  sample. 


1/ 
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f«co i - 1 $•(-;->  j 


p(x,J5)  - J 


r ' '2  2 v • n 


0 


<y3 


■ V+l 


-13  - J 

0 


r10-!  - 


1 SL(Sdk\\ 

2 V 1 fi  ;J 


<yu 


v+1 


-1] 


ir  a-*/ V- 


where  OC  and  V are  constants  for  fixed  n,  which  will  be  explained  in 

a subsequent  paragraph  of  this  section. 

No  difficulty  is  encountered  in  connection  with  0C  curves  (which 

were  obtained  by  numerical  integration)  for  the  acceptance  procedure 

p(x,R)  < p*  °ince  Pr£p  < - Pr-^x  ♦ kS  < U ( whenever  k 

The  procedure  to  accept  if  x + k&  <i  U is  of  course  the  one  commonly- 

used  in  place  of  the  sample  standard  deviation  in  variables  sampling. 

A heuristic  justification  for  the  statistic  p(x,K)  as  an  estimate 

of  p is  the  following.  It  has  been  verified  by  numerical  investigation 

U~-X 

that  the  statistic  ot-s—  is  approximately  distributed  as  non-central  t 

l ^ 

with  degrees  of  freedom  y*  and  eccentricity  7 V +1  whenever and 


V-* 


O'  are  such  that  ~~~  is  sufficiently  large,  i.e.,  when  the  population 
fraction  defective  is  3mall.  The  numbers  oc  and  *•'  are  constants  for 
each  fixed  sample  size  n and  are  determined  as  follows,  under  the 


U-x 


assumption  that  OC  — is  distributed  as  non-central  t^  V~iKp 


K for 


small  p then 


k , „ e 

* 5T  ‘ 


V+l 


H+i 


where  xltl  and  s0  are  the  sample  mean  and  sample  standard  deviation 


/ • j- 


based  on  ^ observations.  Equating  the  first  two  moments  of  the 


— k — 

statistics  x * E and  x^+^_ 


Sy  and  solving  for  oc.  and  v* 


V V+l 

gives  cc  as  a function  of  y1  and  V as  a function  of  n and  k.  In 
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•T>r 

fact  cc  « where  d„  is  defined  by  ER  = d9c7">  and  c9  is  defined  by 

^ C ^ ^ *• 

In  order  to  make  the  constants  independent  of  k the 
values  for  the  limiting  case,  i.e,.  when  p goes  to  zero,  were  chosenr^ 

This  approximation  was  verified  by  the  authors  by  numerical 

U-x 

integration  for  the  statistic  and  reference  to  tables  of  the  non- 
central  t.  It  was  found  to  be  very  good  for  both  large  and  small  n for 
most  of  the  AQL  values  used  in  this  report.  It  is  especially  good  for 
small  n regardless  of  p.  For  the  larger  values  of  p where  n is  large 
the  approximation  did  not  hold  with  great  accuracy.  Hcwever,  it 
could  haj«  been  improved  by  making  and  hence  oLt  functions  of  p as 
well  as  of  n.  This  was  not  done  because  this  would  have  entailed 
increasing  the  number  of  tables  of  estimates,  thereby  complicating  use 
of  the  tables  by  sheer  bulk.  Furthermore,  the  effect  of  this  on  the 
estimate  pXx,R)  would  have  been  slight  since  on  the  average  the  sample 
values  of  fall  in  the  region  of  the  tables  where  the  estimate  p 
changes  very  little  as  n changes. 

Essentially  the  above  argument  indicates  that  p(x,E)  is  approximately 
distributed  as  ^(x^^,s^).  Hence  has  sampling  properties  similar 

to  that  of  p(xy,s^)  where  the  effective  number  of  observations  is  >>  +1. 


This  value  of  obtained  in  this  way  is  exactly  that  obtained  by 

U— x 

Fatnaik  who  suggested  that  ■Trf’  c -s™  may  be  approximated  by  non- 

R 

central  t with  degrees  of  freedom  y*  but  with  eccentricity  fn  -^r. 
This  approximation  is  excellent  and  will  hold  for  all  p but 
would  not  fit  into  the  framework  of  the  proposed  estimate  ^(x,R). 
See  P.  E.  Fatnaik,  The  use  of  mean  range  as  an  estimator  of 

variance  in  statistical  teat,  Biometrika,  vol.  37,  June  1950. 
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